On the Origin of Traveling Pulses in Bistable Systems 
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The interaction between a pair of Bloch fronts forming a traveling domain in a bistable medium 
is studied. A parameter range beyond the nonequilibrium Ising-Bloch bifurcation is found where 
traveling domains collapse. Only beyond a second threshold the repulsive front interactions become 
strong enough to balance attractive interactions and asymmetries in front speeds, and form stable 
traveling pulses. The analysis is carried out for the forced complex Ginzburg-Landau equation. 
Similar qualitative behavior is found in the bistable FitzHugh-Nagumo model. 

PACS number(s): 



Traveling waves far from equilibrium are often formed 
when a uniform state is destabilized by a Hopf bifur- 
cation occurring at a finite wavenumber JlJ. Traveling 
waves or pulses also form from parity breaking bifurca- 
tions of stationary patterns A related mechanism 
that has not received adequate attention involves a par- 
ity breaking front bifurcation in which a stationary front 
solution loses stability to a pair of counter-propagating 
front solutions 0-0]. This bifurcation, sometimes re- 
ferred to as the nonequilibrium Ising-Bloch (NIB) bifur- 
cation, has been found in chemical reactions |7|,|| and 
in liquid crystals [p|JTc|]. Bistable systems, which do 
not necessarily support stationary patterns, may exhibit 
traveling pulses and waves beyond the NIB bifurcation. 
Activator-inhibitor systems with non-diffusing inhibitors 
provide a good example. For fast inhibitor kinetics ini- 
tial domain patterns always coarse grain and converge 
to a uniform state. For sufficiently slow kinetics, and 
beyond the NIB bifurcation, traveling pulses, periodic 
wavetrains, and spiral waves appear. 

Numerical studies of systems with a NIB bifurcation 
indicate that traveling pulses do not appear immediately 
at the front bifurcation point. Instead, there is an in- 
termediate parameter range where initial domains may 
travel but eventually collapse. Only past a second thresh- 
old parameter value do initial domains converge to stable 
traveling pulses ||. In this paper we study the interac- 
tions between a pair of traveling fronts in this intermedi- 
ate parameter range. We find that the balance of repul- 
sive front interactions with attractive interactions and an 
asymmetry between leading and trailing fronts gives this 
threshold parameter value. 



We choose to study the parametrically forced complex 
Ginzburg-Landau (CGL) equation 

A t = (fi + %v)A + (1 + ic!)A xx - (1 + ic 3 )\A\ 2 A (1) 
+ j A* + a , 

where A(x, t) is a complex field and v, c\ , C3 and 7 are real 
parameters. The parameter a can be a complex number, 
but since the final results we present here do not depend 
on its imaginary part we assume a is also real JTT| . The 
CGL equation (7 = a = 0) is often obtained as an enve- 
lope equation for an extended system undergoing a Hopf 
bifurcation at zero wavenumber [[l2j. Then, the variable 
A(x, t) describes weak modulations of the homogeneous 
oscillations. The terms a and •yA* in (Q) represent, re- 
spectively, the effect of parametric forcing with equal and 
twice the system's natural oscillation frequency 
Equation has been introduced recently in the context 
of liquid crystals |l5j . 

The parametric forcing term "/A* breaks the phase 
shift symmetry, A — > Ae 1 ^ , of Eqn. ([!]) and reduces the 
one-parameter family of cw solutions of the CGL equa- 
tion, A = ^ oe i ("-c3M)t+^ ; < <j) < 2tt, to two pairs 
of stable-unstable solutions with fixed <j> values, arising 
in saddle-node bifurcations. Equation (|l|) therefore de- 
scribes a bistable extended system of two stable uniform 
states that oscillate with different phases. The second 
forcing term, a, breaks the parity symmetry (A — > — A) 
of these two states. The front solutions we will be con- 
cerned with connect these two states at x — > ±00. 

A simpler, gradient version of Eqn. (Q) is obtained by 
omitting the linear and nonlinear dispersion terms: 
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FIG. 1. The nonequilibrium Ising-Bloch (NIB) bifurcation 
for front solutions of Eqn. (jl|). For 7 > 7 C there is a single 
stable Ising front with zero speed (solid line). For 7 < 7 C the 
Ising front is unstable (dashed line) and a pair of stable coun- 
terpropagating Bloch fronts appears (solid lines). Parameters: 
fi = 1.0, v = 0.01, ci = c 3 = a = 0.0. 



A t = nA + A xx - \A\ 2 A + 7A* + a . 



(2) 



A physical application of (g) is Rayleigh-Benard con- 
vection with periodic spatial modulation of the cell 
height jl6| or heating jL7|. When a — 0, Eqn. (g) has 
three types of stationary front solutions 

I(x; a) = crA tanh (—=A Q x) , (3) 
V 2 

B±(x] a) — <jAq tanh(/cx) ± iy/ji — 37 sech(fca:) , (4) 



where Ac = \ffjb + 7, k = */2rj and a = ±1 is the 
front polarity which stems from the reflection symmetry 
x — > —x of Eqns. (||) and ([!]). The front solutions I(x; a) 
and B±(x;a) are equivalent to Neel (Ising) and Bloch 
domain walls in ferromagnets with weak anisotropy ^] 
and will be referred to as Ising and Bloch fronts. The 
Ising front I(x; a) loses stability as 7 is decreased past 
the critical value 7 C = /j/3. At that point the two Bloch 
fronts B±(x; a) appear and are stable (J]. 

The non-gradient terms associated with is, c\ and C3 re- 
move the degeneracy of the three stationary solutions (||) 
and ([IJ). With any of these terms nonzero, the two Bloch 
fronts propagate in opposite directions at a speed pro- 
portional to the corresponding coefficient, 2/, c% or C3 [|J. 
In that plot of the front velocity c, versus 7 yields 

the NIB bifurcation diagram shown in Fig. [|. 

To study front interactions we consider the symmetric 
(a = 0) and nearly gradient case, where front solutions 
of (Q) can be expanded around front solutions of the gra- 
dient system (||) . We introduce a small parameter e<l 
and assume that the constants v, a, c\, and C3 are all of 
order e. We also assume proximity to the Ising-Bloch 
bifurcation point, /i — 37 ~ yfi. A traveling domain so- 
lution of Eqn. (pi) is sought as 



A(x,t) =B + [x-xi{T);+l]+B_[x-x r (T);-l] (5) 
-A + R{x, T) , 

where x r and xi are the positions of the leading (right) 
and trailing (left) Bloch fronts, T = et is a slow time, B± 
are given by (Q) , and R is a small correction term of order 
e. The two polarities (a ~ ±1) are necessary to construct 
a domain bounded by the fronts. The two types of Bloch 
fronts, _B_ and B + , make the domain traveling instead 
of shrinking or expanding. We assume that the domain 
is much wider than the width of the fronts, k , or more 
accurately, that exp[— 2k(x r — x{)) ~e<l. 

Our objective is to derive equations of motion for the 
front positions, Xi and x r . Using the anstaz (^|) in (|l|) and 
following the method of Refs. ]lq] , we find the following 
solvability condition for the right (leading) front: 
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where c.c. stands for the complex conjugate. A simi- 
lar condition is obtained for the left (trailing) front. In 
Eqn. (^|) the dot over x r represents the derivative with re- 
spect to the fast time t, and Af is a nonlinear interaction 
term 

Af=(B*_+ B* + )(A - B+){B_ - A ) 
+ (B_ + B + )(A -B + )(B*_-Ao) 
+ {B- + B + )(A - B* + ){B- - A Q ) , 

where the arguments of B± are as in (||). Analytical eval- 
uation of the nonlinear interaction integral in (ra ) leads 
to the following equations, with r\ — y/ji — 37 Pq j: 
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2a + ixr\[-v + c 3 ^ + (c a - 03)7] 

- 8A 3 e- 2k{x — Xl) + iAorfe-^"-^ 



(7) 



-kA x t = ~2a + irr][-v + c 3 /x + (ci - €3)7] 

+ 8A 3 e- 2k{Xr - Xl) - iAo^e-^"-^ . 

Combining Eqns. gives a single equation of motion 
for the distance between the two fronts, L 



kA L = 3a- 12Ale 
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The first term on the right hand side describes the effect 
of the broken symmetry between the two Bloch fronts; 
the initial domain expands (shrinks) in time when the 
leading front is faster (slower) than the trailing front. 
The second term describes an attractive front interaction 
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generated by the real parts of the Bloch front solutions. 
The last term, generated by the imaginary parts of the 
Bloch front solutions, describes a longer range repulsive 
interaction. The repulsive interaction strengthens as 7 is 
decreased below the Ising-Bloch bifurcation point, 7 C = 
/Lt/3, and becomes dominant at sufficiently small 7 values. 

Notice that in Eqn. (||) the non-gradient terms asso- 
ciated with is, ci, and C3 have disappeared. To leading 
order the effect of these terms is to grant the two Bloch 
fronts equal speeds. Therefore the distance between 
the fronts is not affected. This suggests that the exis- 
tence and stability of pulse solutions (moving domains) 
of Eqn. (0) can be studied starting from the gradient ver- 
sion (||). The latter can be written as A t = -8H/5A*, 
where S/SA* stands for the variational derivative, and 
the Lyapunov function (pseudo-Hamiltonian) is 

/+00 
Hdx , (9) 
-OO 

H = \A X \ 2 -»\A\ 2 + ±\A\* 

-l 7 [A 2 + (A*) 2 ]-a(A + A*). 

This gradient representation implies relaxational dynam- 
ics toward a minimum of H . The part of H which de- 
pends upon the separation distance L between the Bloch 
fronts can be found following Ref. |jl9). The resulting 
effective pseudopotential of the interaction for /1 ps 37 is 

(10) 

The extrema of U (L) coincide with stationary solutions 

To find traveling pulse solutions of O), which corre- 
spond to stationary (fixed point) solutions of (||), we set 
L = in (|) or dU/dL = in @ and solve the result- 
ing quadratic equation for z = exjp(-kL). The solutions 
are: 

L = -AT 1 In (?? 2 ± V^ 4 + 4A)a) + In 2A . (11) 

In the symmetric case, a = 0, and for rj > 0, there is 
only one (finite) pulse solution. A linear stability anal- 
ysis indicates the solution is unstable. This result com- 
plies with an earlier finding reported in Ref. |2(i(| . For 
a > the leading front is faster than the trailing front. 
Again, only one pulse solution exists and a linear sta- 
bility analysis indicates the solution is unstable. The 
same conclusion follows from a graph of Eq. (|l0|): the 
single pulse solution corresponds to a maximum of the 
interaction pseudopotential U(L). For a < no pulse 
solutions exist if 77 4 < 4Ao\a\. A saddle- node pair of 
pulse solutions appears at 7 = 7 P (a), where 7p(a) solves 
(/j — 37 P ) 2 = 4^//j + 7 P |a|. Graphs of these solutions in 
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FIG. 2. The distance, L, between the front and back of 
a pulse solution for 7 near the NIB bifurcation. The solid 
and dashed lines represent the stable and unstable branches 
solutions from Eqn. (0). The crosses are data from direct 
numerical solution of Eqn. (jl|). Parameters: y, = 1, v = 0.01, 
a = -0.001, ci = c 3 = 0. 

the L — 7 plane are shown in Fig. |^. The upper and 
lower branches represent stable and unstable solutions, 
and pertain to a minimum and a maximum of U (L) , re- 
spectively. Also shown in Fig. || are results from direct 
numerical solutions of Eqn. p) showing the stable pulse 
branch. The agreement for small 77 is within 5%. The 
shape of the stable traveling pulse is shown in Fig. |^a. 

The conclusion that no stable pulses exist for a > 
is a result of the specific ansatz (||) for an "up" pulse as 
shown Fig. ||a. From the symmetry A — > —A of (|l|), a 
symmetric stable "down" pulse exists for a > 0. The 
shape of this pulse is displayed in Fig. ||b. A phase di- 
agram of all stable front and pulse solutions is shown in 
Fig. [§i. Three main regions can be distinguished: (I) 
The entire Ising regime, 7 > 7 C = ^(/3, where only a 
stable Ising front solution exists. In this region domains 
shrink or expand but do not travel or form pulses. (II) A 
region in the Bloch regime, r y p {a) < 7 < j c , where a pair 
of stable counterpropagating Bloch fronts exist. In this 
region domains travel but do not form pulses; they cither 
expand to infinity or shrink and collapse. (Ill) The rest 
of the Bloch regime, 7 < 7 P (a), excluding the a = line 
(the 7 axis), where stable traveling pulses exist in addi- 
tion to the pair of Bloch fronts. The 7 axis separates 
regions of up and down pulses. 

This behavior is rather general. Shown in Fig. ^b is 
an analogous phase diagram for the bistable FitzHugh- 
Nagumo (FHN) model 

u t — u — u 3 — V + u xx , 

v t = e(it - a±v - a ) + Sv xx , (12) 

where u and v are real fields and e, 5, a\ and ao are real 
constants. The parameters e and ao play a similar role 
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FIG. 3. Pulse solutions to Eqn. (|lj) near the NIB bifurca- 
tion, (a) A stable "up" pulse, a < 0. (b) A stable "down" 
pulse a > 0. 



as 7 and a, respectively, and the same three regions (I), 
(II) and (III) appear in the e-ao plane. 

The existence of an intermediate parameter region 
where traveling domains collapse rather than converge 
to stable pulses, has been observed in numerical sim- 
ulations of model equations p[,pT[. It is also ex- 
pected to be observed in a number of experimental 
systems including catalytic reactions on platinum sur- 
faces pi | 
fields 



liquid crystals subjected to rotating magnetic 
|, oscillatory chemical reactions subjected 
to periodic forcing ^3|, and the Ferrocyanide-Iodate- 
Sulfite reaction |E4|. 
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